Observation of Geometric Phases for Three-Level Systems using NMR Interferometry 
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Geometric phase (GP) independent of energy and time rely only on the geometry of state space. It 
has been argued to have potential fault tolerance and plays an important role in quantum information 
and quantum computation. We present the first experiment for producing and measuring an Abelian 
geometric phase shift in a three-level system by using NMR interferometry. In contrast to existing 
experiments, based on the geometry of S 2 , our experiment concerns the geometric phase with the 
geometry of SU(3) /U(2). Two interacting qubits have been used to provide such a three-dimensional 
Hilbert space. 
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When a quantum mechanical system evolves cyclically 
in time so that it returns to its initial physical state, its 
wave function can acquire a geometric phase factor in 
addition to the familiar dynamic phase [1]. If the cyclic 
change of the system is adiabatic, this additional factor 
is known as Berry's phase [2]. Otherwise, it is related to 
Aharonov-Anandan (AA) phase [3] that has been pointed 
out to be a continuous version of earlier Pancharatnam 
phase [4]. 

Geometric phases (GP) independent of energy and 
time rely only on the geometry of state space. It is there- 
fore resilient to certain types of errors and suggests the 
possibility of an intrinsically fault-tolerant way of per- 
forming quantum gate operations [5, 6, 7]. This potential 
value makes it important to observe and further apply 
GP in different quantum physical systems. The observa- 
tions of GP began from earlier spin-polarized neutrons 
through a solenoid [8] , polarized light through a helically 
twisted optical fibre [9] , and a pair of coupled protons in 
magnetic field using NMR [10] to the recent supercon- 
ducting qubit experiment [11]. The principle of them are 
usually the same. That is, in a two-level state space the 
geometry of it corresponds to a sphere S 2 and the GP 
equals to one half the solid angle subtended by closed 
paths on S 2 . 

When one generalizes to a three-level quantum sys- 
tem [13, 14], the geometry of S 2 gets replaced by a 
four-dimensional geometric space SU(S)/U(2) or part of 
sphere S* 7 . Then evolutions of state correspond to actions 
of SU(3) on SU(3)/U(2) that is different from that of 
SU(2) on S 2 for the two-level case. In order to observe 
GP, one way to vanish the dynamical phase is closely 
linked to the geodesic in ray space (see below). For two- 
level case, the geodesic in ray space happens to coincide 
with that on S 2 . In contrast, it is a plane curve instead of 
geodesic on S 7 for three-level case. The GP for any cyclic 
evolution in three-level ray space are no longer related to 
solid angles on S 7 but referred to Bargmann invariants 
[15, 16]. All of these differences indicate the observation 



of three-level GP technically more difficult [18]. 

In this letter, we report an experimental observation 
of three-level GP by using NMR interferometry. The 
three levels referred in the experiment are chosen from 
a two spin- 1/2 interacting system. Unitary evolutions 
for implementing cyclic paths in the three-dimensional 
ray space are ensured by quantum controlled logic gate 
operations [12]. Aimed at obtaining a measurable GP, we 
evolve the target state while keeping the reference state 
unchanged to produce a relative phase between them. 
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FIG. 1: An illustration of the parameter space for all three- 
level states. The local coordinates 9, <f>, xi and X2 are such 
that (6, 4>) define a point in the positive octant of S 2 and 
(Xi)X2) for (0, <fi) fixed define a point on a torus. 

In the three-dimensional Hilbert space 7i 3 , an arbi- 
trary state can be expressed as 
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where the real parameters have the range 9,(f> £ [0, ?] 
and xi>X2 S [0, 27r). It has an one-to-one correspon- 
dence, omitting a global phase 77, to the point on an oc- 
tant of S 2 plus a torus (Fig. 1). Consider a state evolves 
from \tp(si)) to \ip(s2)) with si,2 being the curve param- 
eters determined by the Eq. (1). Corresponding to this 
evolution, in H? there is a continuous piecewise smooth 
parametrized curve, C — {^(s)|si < s < S2}, and its im- 
age in the ray space 1Z is likewise continuous and piece- 
wise smooth denoted by C = {p(s) — \ip(s)){ip(s)\\si < 
s < S2}. Then the GP (3 associated with the cure C 
equals to the difference between a total phase ftot and a 
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dynamical phase jd [14], that is, 

0[C] = ^tot[C]- 7 d[C], 
Vtot[C] = arg('0(si)|V'(s2)), 



(2) 



7d[C] = - 



d 

ds(i/j(s)\i—\ip(s)), 



with both (y9 t ot and 7^ being functionals of the curve C. 
If the curve C is closed, the state change can be simply 
expressed as \tp(s2)) = exp{i(7d[C] + /3[C])}\tp(si)). 

The geodesies in ray space 1Z are given through vari- 
ations of a nondegenerate positive definite length func- 
tional [16]. In two- level systems geodesies are related 
to the parellel transport condition. For the three-level 
case, every geodesic in ray space has a vanishing geomet- 
ric phase and it plays a crucial role in the observation 
of geometric phases in the following. The simplest de- 
scription of geodesic can always be achieved as follows 
[14]. Let pk and Pk+i denote the end points of a smooth 
curve C associated with unit vectors ipk and tpk+i in W 3 . 
There is a requirement for the chosen state vectors that 
(ipkl'i'k+i) must be real positive. Then the geodesic C geo 
connecting pk to pk+i is the ray space image of the curve 
Cgco = {"0( s fc)|O < s k < s°} and 



1p(Sk) = i>k COSSfc + 



tpk+i - tpki^kl^k+i) 

V 1 ~ {^k\lpk+l) 2 



sinsfc, (3) 



with < Sk < s° and s£ = a,rccos(ipk+i\ipk) ■ From the 
Eq. (3), one can see that ^j(O) = -0fe and |i/'(s°)}(V'(sfc)| = 
\ip k+ i)(ipk+i\ = Pk+i- For a set of points pi,p 2 , • ' ' , Pn C 
1Z in order, suppose that no two consecutive points 
are mutually orthogonal and that p n and p\ are also 
nonorthogonal. So we can obtain a closed curve C in 
1Z in the form of a n-sided polygon by joining these n 
points cyclically with geodesic arcs. The geometric phase 
is then according to the Eq. (2) 



f3[C] = arg(Vi|^)-arg(Vi|^ 2 ) 
= -argTr(pip 2 • • ■ Pn), 



- arg(^ n |V'i) 



(4) 



in which it has used relations of = = Pi 

and p\—p\. The Eq. (4) combined with geodesic con- 
dition, i.e., (ipk\ipk+i) is real positive, shows a vanish- 
ing dynamical phase for these cyclic evolutions. It thus 
provides us a convenient evolution way to observe the 
geometric phase. 

Experiments are performed on the three-dimensional 
subspace of two interacting spin-i nuclei — spin a ( 13 H) 
and spin b ( 1 C) in the 13 C-labeled chloroform molecule 
CHCI3. The reduced Hamiltonian for this two spin sys- 
tem is, to an excellent approximation, given by H = 
+ u b I b z + 2-kJIII\. The first two terms in the 
Hamiltonian decribe the free precession of spin a and 
spin b around the magnetic field Bo with Larmour fre- 
quencies w q /2tt w 400 MHz and w b /27r w 100 MHz. 
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FIG. 2: Energy level diagram for (solid lines) two spins cou- 
pled by a Hamiltonian of the form of 2nhJlll^ and (dashed 
lines) two uncoupled spins. 



The third term of the Hamiltonian describes a scalar 
spin-spin coupling of the two spins with J = 214.5 Hz. 
Experiments were performed at room temperature on a 
Bruker AV-400 spectrometer. If we denote the spin up 
and down by |0) and |1), the energy levels of such sys- 
tem are displayed in Fig. 2. It has four levels written as 
{|00}, 1 01 ) , 1 10) , 1 11)} corresponding to energy eigenvalues 
{^h(—u>i - L02 + nJ), \h{-wi + 0J2 - 7rJ), - a>2 — 

7rJ), \h{u)\ + uo 2 + nJ)}- We choose basis states { 1 00) , 
1 10) , |H)} to construct the desired three-level space Tt 3 
and 1 01) as the reference state which keeps unchanged 
during evolutions. 

The system was first prepared in a pseudopure state 
|00) using the method of spatial averaging [17] with the 
pulse sequence 



^(tt/3) - G z - J&7T/4) -> — - J$(tt/4) 



G, 



(5) 

which is read from left to right (as the following se- 
quences). The rotations (angle) are implemented 
by radio-frequency pulses. G z is a pulsed field gradient 
which destroys all coherences (x and y magnetizations) 
and retains longitudinal magnetization (z magnetization 
component) only, represents a free precession period 
of the specified duration under the coupling Hamiltonian 
(no resonance offsets). 

The complete sequence started by preparing the ini- 
tial superposition state ^(|00) + 1 01) ) with a Hadamard 
operation on the second qubit of the pseudopure state 
|00). Then the reference term |01) was kept un- 
changed through bipartite control operations as shown 
in Fig. 3. The |00) term (denoted by 
was first evolved to ^2) = coss^OO) + sin 1 10) 
with unitary operation Uf(si), then to state ^3) = 
(cos S® cos — e l6 sin sin s° c °s <p) |00) + (sin cos + 
e cos S] 1 sin s° cosy) 1 10) + sin^sins^lll) with f/f(s2), 
and last to state \ipi) = e l/3 |-0i) with f7f(s3). Corre- 
sponding to three smooth geodesies, the unitary opera- 
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tions can be factored into more clear form 



where 



= R(si), 

= i?(s°)i? 23 (6>,0) R(s 2 ) R^ 1 (9,tp,0)R- 1 (s° 1 ), 



(6) 



R{s k ) 



( cos Sfc — sin Sfc 0\ 
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and the 5J7(2)23 subgroup element 
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Parameters £,Xj t j s 3 in Eq. (6) are determined 
by the reparametrization ^3) = e** cos S3 1 00) + 
gHJ+x) sin s§ cost) 10} + sins§ sinrjll) and curve param- 
eters Sfc (k = 1, 2, 3) have ranges < s^ < s" . Obviously 
the chosen unit vectors ipk and tpk+i satisfy the condition 
of (ipk\ipk+i) being real positive. This combined with the 
Eq. (4) shows a vanishing dynamical phase during these 
cyclic evolutions and we obtain the GP 



(3[C] = arg(coss5 cossj — e z9 sin s? sin S2 cosy). 



(7) 



So after one cyclic evolution described above, it effec- 
tively produces a GP and can be measured as a relative 
phase shift between |0)& and 1 1 ) t, for the qubit b, i.e., 
^(|00) + |01» - ^(etf|00) + |01)) - |0) 8 ®^(e^|0) + 
\l))b- At last the local phase f3 can be read out directly 
by a phase sensitive detector on qubit b in NMR. 
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FIG. 3: Experimental network: two spin-1/2 nuclei perform 
unitary evolutions controlled by each other. Each circle at the 
second line means that performs its linked unitary evolution 
when the second nucleus at |0) state. Each dot at the first 
line means that performs its linked unitary evolution when 
the first nucleus at |1) state. 

In Fig. 4 we show the measured phase and its de- 
pendence on different points A, B, C in ray space char- 
acterized by parameters {s?, s°, 0, tp}, all carried out at 
6 = 7r/4, and total pulse sequence time T for cyclic evo- 
lution is about 5 ~ 25ms for different evolution path. In 
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4: Experimental results on the GP f3 versus the paramo- 
Changing si, s® means changing the positions of 
the points B and C. The evolution paths have been depicted 
out in parameter space and the theoretical curves are marked 
out by lines, (a) It shows the result in the case of 6 = 7r/4 
and tp — 0; the ABC on the octant of S 2 is a curve while it 
runs a period on the torus, (b) It shows the result in the case 
of 9 — 7r/4 and tp = 7r/4; the ABC on the octant of S 2 is a 
triangle while it runs a period on the torus. 



Fig. 4 (a) and (b), we set tp — and tp = n/4 respec- 
tively with which geodesies have disparate trajectories in 
ray space. The measured phase is in all cases seen to 
fit the theoretical curve (7) well with a root-mean-square 
deviation across all data sets of 5.6 degree. Thus, all re- 
sults are in close agreement with the predicted geometric 
phase, and it is clear that we are able to accurately con- 
trol the amount of phase geometrically. 

The controlled operation between the two qubits plays 
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the main role in the experiment. It goes as that the qubit 
a (or b) undergoes a SU{2) operation if the qubit b (or 
a) is in |1) while kept unchanged if it in |0). This is 
used to realize the controlled operations R and i? 2 3 • The 
concrete operations go as follows. 

For the subsystem of qubit a, we can write the reduced 
Hamiltonian 

H a = ujj a z + 2irJm» z I a z = [w a - 2nJ(d b - \)]I a z , 

where m b z is the eigenvalue of I z (= ±i) and d b the cor- 
responding computational value (= 0, 1). If we use a 
rotating frame with a frequency of ui' a — uj a and u)' b = Wb, 
the Hamiltonian turns into, for d b = 0, = tt,JI z , 

while it becomes = —irJI z for d b = 1. This Hamil- 
tonian generates controlled rotations around the z-axis. 
To generate the control gate R(Sk), we rotate the rota- 
tion axis using radio-frequency pulses. To generate a 2Sk 
rotation around the y-axis, e.g., we use the sequence 

This represents the controlled gate operation R(Sk). 

For another controlled gate operation i?23(x> T i ~ £)> we 
have to reverse the roles of control and target qubit and 
apply the following sequence to qubit b: 

Rl(-<P) - RK-tt - 0) - ^ -> Rl(n + 13) 
^R 2 M^Rl(a,f3,<j>), 

Rn(a,/3,(f>) denote to rotate the second qubit a around 
the axis n(/3, (/>), and a, /?, <f> is calculated from \, r, — £. 

In conclusion, when a quantum mechanical system 
evolves cyclically in time so that it returns to its initial 
physical state, its wave function can acquire a geometric 
phase factor in addition to the familiar dynamic phase. 
Geometric phases (GP) independent of energy and time 
rely only on the geometry of state space. It is therefore 
resilient to certain types of errors and suggests the possi- 
bility of an intrinsically fault-tolerant way of performing 
quantum gate operations, we present the first experi- 
ment for producing and measuring an Abelian geometric 
phase shift in a three-level system by using NMR inter- 
feromctry. In contrast to existing experiments, based 
on the geometry of S 2 , our experiment concerns the ge- 
ometric phase with the geometry of SU(3)/U(2). Two 
interacting qubits have been used to provide such a three- 
dimensional Hilbert space. 
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